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The behaviour of ferromagnetic systems with single-ion anisotropies and magnetic fields in more 
than one direction is investigated with many-body Green's function theory generalizing earlier work 
with uniaxial anisotropies only. It turns out to be of advantage to construct Green's functions in 
terms of the spin operators , and , instead of the commonly used S'^ , S~ and 5*^ operators. 
The exchange energy terms are decoupled by RPA and the single-ion anisotropy terms by a general- 
ization of the Anderson-Callen decoupling. We stress that in the derivation of the formalism none of 
the three spatial axes is special, so that one is always able to select a reference direction along which 
a magnetization component is not zero. Analytical expressions are obtained for all three components 
of the magnetization and the expectation values ((5*^)^), {{S")'^) and {{S^)^) for any spin quantum 
number 5*. The formalism considers both in-plane and out-of-plane anisotropies. Numerical calcu- 
lations illustrate the behaviour of the magnetization for 3-dimensional and 2-dimensional systems 
for various parameters. In the 2-dimensional (monolayer) case, the magnetic dipole-dipole coupling 
is included, and a comparison is made between in-plane and out-of-plane anisotropies. 



I. INTRODUCTION 



The spontaneous magnetization of bulk ferromagnets has been successfully described by a Heisenberg model within 
the framework of many-body Green's function theorji^, where the Tyablikov decoupling (random phase approx- 
imation: RPA) is the crucial approximation for the exchange energy term. Usually, only the z-component of the 
magnetization (S^) was considered. Calleni^ developed a formalism which lead to a closed-form expression for the 
magnetization (S^) in ferromagnetic (FM) systems for any spin quantum number S. This theory has been helpful in 
studying various systems described by Heisenberg HamiltoniangiiSiSii. 

Recently, experimental work has provided an incentive to calculate more than one component of the magnetization. 
For example, under the influence of an external field perpendicular to the anisotropy axis, the magnetization rotates 
with the variation of the field strength. When a FM layer is covered by an anti-ferromagnetic layer, the magnetizations 
across the interface may have different orientations^. Moreover, the antiferromagnetic sublattice magnetization may 
be non-collinear— . In FM ultrathin films, the magnetization may rotate when the temperature or film thickness are 
varied; see e.gii2iii*i2ii^. 

Frobrich et ali^ii^ have developed a method using Green's function theory to calculate more than one component 
of the magnetization for FM films. They investigated the reorientation of the magnetization of FM films caused by an 
external field and the dipole-dipole coupling by applying RPA to the exchange energy term and the Anderson-Callen 
decouplings^ to an out-of-plane single- ion anisotropy. Later an exact treatment of the single- ion anisotropy term 
was proposed by going to higher-order Green's functions^^. This exact treatment and a comparison with Quantum 
Monte Carlo calculationsi^ showed that the Anderson-Callen decoupling is a good approximation when the single-ion 
anisotropy is small compared to the exchange interaction. The Anderson-Callen decoupling can therefore often be 
a useful approximation in this situation, particularly as it is much simpler to apply than the exact treatment. In 
Refiii, general analytical formulas were derived on the basis of RPA and the Anderson-Callen decoupling for the three 
components of the magnetization for 3D and 2D FM systems for any spin quantum number. The Green's function 
formalism was also applied to uniaxial out-of-plane and uniaxial in-plane exchange anisotropiesiS*i2iS£. 

The papers above [12-17] dealt only with a uniaxial out-of-plane single-ion anisotropy. In the present paper, we 
generalize the model by introducing single-ion anisotropies in all directions of space, thus obtaining the in-plane and 
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out-of-plane anisotropies as special cases, and also allowing for a reorientation of the magnetization from out-of-plane 
to in-plane and vice versa, depending on the chosen parameters. The new key features are the formulation in terms of 
Green's functions defined in the spin operators 3^,8^, instead of the ,S~ operators, and the generalization 
of the Anderson-Callen decoupling to the terms coming from the single-ion anisotropies for all directions in space. 
The new formulation is advantageous if several non- vanishing components of the magnetization are present. 

The paper is organized as follows. Section 2 presents the general formalism for any spin quantum number. Section 
3 shows numerical results for the 3D case and for a monolayer film (2D). For the latter case we also discuss the effects 
of the magnetic dipole-dipole coupling. Section 4 contains our conclusions. 

II. THE FORMALISM 

The Hamiltonian for a ferromagnetic system is 

= E • - E[^2.(5f )2 + K^yiSjf + K^ASm - E B ■ S. . (1) 

The first term is the Heisenberg exchange interaction (J > 0). The second term describes single- ion anisotropies in 
all three directions of space where we assume that the anisotropy constants are greater or equal to zero, > for 
a = x,y, z. Note, that due to the exact relation (S^)'^ + {S^Y + ('S'^)^ = S'(S'-|- 1) only two of the are independent. 
The third term is the Zeeman energy, arising from an external magnetic field. 

We use retarded many-body Green's functions defined in accordance with Bogolyubov and Tyabliko\s2i as 

G^At " i') = {{A,:B,))^ = ~i6{t ~ mA,B,],), (2) 

where the operators Ai and Bj depend on the problem under consideration, and the subscripts i and j label lattice 
sites. 9{t — t') is the Heavyside function, and (...) denotes the thermodynamical expectation value. Furthermore 

[A,, Bj]„ = A,B.j + T] BjA, = i^y- „, r] = ±1. (3) 
After a Fourier transformation to energy space, we obtain the equation of motion: 

cj{(A,- B,))r,,^ = {[A^.B,]r,) + {({A, H]; B,)),,,^. (4) 
The Green's function is further Fourier transformed to momentum space 

ij 

The aim is to obtain the Green's function g{]<.,uj) from a solution of the equations of motion. Then, statistical 
averages of the product of the two operators Bj and Ai (correlation functions), may be cast in terms of correlations 
in momentum space, Ck: 

{B,A,)^^J2^'''''^'''^^^- (6) 

k 

The Ck can be calculated by the well-known spectral theoremi 

^^ = V- I ^vfc^^[.9(k,c.-^z0+)-.g(k,..-z0+)], (7) 
27r J e^i'^^-^ + 77 

where T is the temperature and fc^ is the Boltzmann constant. 

In order to solve the equations of motion, the higher-order Green's functions appearing in Eq. (4) have to be 
decoupled. The following derivations make use of Fermionic Green's functions (77 = -1-1) instead of the Bosonic 
Green's functions (77 = —1) of earlier work, where zero eigenvalues of the equation of motion matrix have to be 
treated carefully, e.g>i2iS. 

In previous work [12-19], the operators A and B were the spin operators: 



A = 5"^, K = +, -, z, and B = (5-)™(5^)", 



(8) 
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where m and n are zero or positive integers, necessary for treating larger spin values S. The z-axis was chosen to be 
the particular direction in which the single-ion anisotropy is active. 

In this paper, in order to to treat the three directions in space on an equal footing, we use the following spin 
operators to construct the Green's functions: 

A = 5", a = X, y, z, and B = (S'^)'(S'^)"(5^)", (9) 

where l,m,n are integers or zero. 

The advantage of such operators is that the S°"s are equivalent to each other. A particular direction in space is 
not preferred. 

Analogous to the usual formulation, wc apply RPA to the terms coming from the exchange interaction in the new 
variables. This corresponds to the following factorization: 

((Sr^f ; B,)) = (5r)((5f ; B,)) + (5f )((5f; B,)), a, p = x, y, z, and * ^ I. (10) 

RPA is not an appropriate decoupling for the on-site terms coming from the single-ion anisotropy. For these, we use a 
generalization of the Anderson-Callen decoupling (originally derived only for the uniaxial out-of-plane anisotropies^) 
by assuming that the decoupling form factor $q, has an identical form in all spatial directions (a = x, y, z). 



{{S^S^ + S^S'^-B,)) = (5f )$„((^f = (5f)$^((5f a, (3 ^ a;,y,z and a^p. (11) 

Here 
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" = 2 ( 1 - 7^[S{S + 1) - {ST St)] ) , a = X, y, z. (12) 



As appropriate anti-commutator Green's functions we choose 



G= I ((5f;i3,))+i I . (13) 



1 



Applying the decoupling procedures (10, 11), we obtain an equation of motion from which the Green's functions can 
be determined 

/([5f;i?,]+i)\ 

(^I-P)g = F^,=. {[Sl,B,U) . (14) 
\{[Sl,B,U)) 

One of the main advantages of the present formulation is that the matrix P turns out to be Hermitian, whereas in 
earlier workiStiii^iii, the corresponding matrix at this stage of the decoupling was real and non-Hermitian. 

' iH^ -iHy \ 
-iH^ iH^ \ . (15) 
iHy ~iH, j 

Here 

-ffa = (5")(Jo- Jk) + S„ + i^2„(^")$„, a = x,y,z. (16) 
The quantities J's are defined as 

Jk = e* '"^ and Jq = 7 J, (17) 

where the summation is over the nearest neighbours and 7 is the number of nearest neighbours. 
Because the matrix P is Hermitian, it has real eigenvalues: 



ioi = 0, ^2,3 = ±JHI + m + i72 ^ ±E^. (18) 
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The eigenvector matrix can be determined analytically, where the columns correspond to the eigenvalues 0, — -Bk, +E]i, 
respectively 

hx/hz —{hz — ihxhy) / {hx + ihyhz) —{hz+ihxhy)/{hx — ihyhz 
U=l hy/hz -i{hl + hi) /{hx + ihyhz) i{hl + hi) /{hx - ihyhz) I, (19) 



and its inverse is 



where we have defined 



^ ; 2hxhz ^hyhz 2hj. 

= - -hxhz ~ ihy -hyhz + ihx hi + h^ I , (20) 
* -hxhz + ihy -hyhz - ihx hi + hy ' 



ha ^ a = x,y,z. (21) 



The Green's function is obtained hyi^ 



=EE ^^^+1 =E i^-^^'-^" . (22) 



A=lT = l ^ A=l 



Applying the spectral theorem (7), we get the correlation functions (with r/ = +1) in momentum space 

/ {BS-) \ 

Ck = (BSy) = RkF+i,k. (23) 
V (BS^) J 

The index k is introduced to indicate that the inhomogeneity for the anti-commutator is momentum-dependent, 
whereas that for the commutator is not. 

Now we exploit the relation between the anti-commutator and commutator inhomogeneities, 

F+i,k = F_i + 2Ck, (24) 

to obtain 

(I-2Rk)Ck = RkF_i. (25) 

This equation is written explicitly as follows 



/ ihz —ihy \ 






-ihz ihx 






\ ihy —ihx 1 






coth{Ei,/2kBT) 


-ihz 


ihy 
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-1,2 




-1,3 



ihz coi\i{E^/2kBT) -ihx 1 I F^i^2 1 . (26) 

—ihy ihx coth(i?k/2fcBT) 

It should be pointed out that these equations are not independent of each other: On the left hand side of equations 
(26), the first row multiplied by hx plus the second row multiplied by hy plus the third row multiplied by hz is equal 
to zero. Applying the same condition to the right hand side, we find 

+ hyF_i^2 + = 0. (27) 

Taking B of Eq. (9) as = 5^, and S*^, respectively, we get 

hy{S') = hz{Sy), hx{S') = hz{S^) and hxiS^) = hy{S^). (28) 

Equations (27) and (28) are called regularity conditions because they are equivalent to similar equations in Refii^, 
which were derived there from the condition that the commutator Green's functions are regular at the origin, w = 0. 
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For the convenience of the discussion below, we now define two quantities Qa and hap as 

1 coth(^k/2fcBT) K Ha . 

^°'^1^t1^ 1 and /la/3 = — = —, a,l3 = x,y,z. (29) 

k ^ 

Before further discussion, we recall that Gallon'^ derived a general expression for calculating (S*^) only: 

_ (5 - Q)(Q + 1)^^+1 + (^ + 1 + Q)Q^^+i 

(Q + 1)2S+1 _ Q2S+1 ' ^''"^ 

where 5 is the spin quantum number and Q — ij {^Q^{^)/kBT _ j^-j -^^jth (jj(k) being the dispersion relation. We 

found^^ that putting Eq. (30) into a symmetric form has some formal advantage 

_ (2^ + 1 - + 1)^^+^ + (25 + 1 + Qi)(Qi - 1)^^^+^ 

2[(Qi + l)2S+i_(g,_i)2S+i] ' ^-^^^ 

where Qi = ^ coth(cj(k)/2fcBr) = 1 + 2Q. 

Now we are prepared to derive the required formulas for the components of the magnetization from Eqs. (26) 
and (27). We deal first with the case of an in-plane magnetization, say in the a;y-plane. It is comparatively easy to 
get general expressions for this case. Then we proceed to the most general case, where it does not matter in which 
direction the magnetization points. 

2.1. Magnetization in the sy-plane. 

Here we take the magnetization to be in the xy-plane. In this case {S^) ~ and = 0, and wc assume that 
hx ^ 0, i.e. we select the a;-direction as reference direction. Two of the Eqs. (26) can be rewritten as 

—2hyxCi + 2C2 = hy^F^ii — F_i2 + iQxF-1,3- (32) 

These equations allow us to obtain the expressions for the magnetization and all correlations for each spin value S 
by specifying the proper B — {S^Y {S^)"^{S^)'^ operators, and using the regularity conditions and some spin algebra. 
Because of the Anderson-Callen decoupling (11, 12), the correlations ((5^)^), ((5*^)^) and ((5^)^) are also needed. 
We list in Table I the explicit expressions for (S^), ((S*^)^) and ((5^)^} for spin quantum numbers S = 1/2 to 5/2, 
where 

P2 = 1 + hi,. (33) 

Note that this quantity is momentum independent due to the regularity conditions (28), hyx = {S'^)/{S^). 
Table I. The analytical expressions (5^), ((S'^)^) and {{Sy)^) and for 5 = 1/2 to 5/2. 

5 (5") ((5^)2) {{sy)^) 

1/2 1/2Q, 1/4 1/4 

4Q^ 2(Ql + l) 2iQl+hlJ 

bQl+P^ + 5Ql+R^ + 12hl^ 

■2QAQl+P^) 4(Q2+P2) 4(Q2+p2) 



1 

3/2 
2 



4Q.(5Q^+3P") WQl+i2Ql+f>P^+&QlP^ WQl+i2Qlhl^ +6P" feg, +6Ql P^ 



5Qt+10Q'iP'-'+P'^ 5Qi + WQ'iP'^+P^ SQ^ + lOQiP^+P" 

r/r, 35Qt:+42QlP^+3P* 35Q^ + 14Q^(3P^ + 16)+3(P"+32)P'' 35Q^+14Q^(3P^ + 16feg^)+3(P^+32feg^)P' 
/ 2QlpQfTT0QjP'T3iF) 4(3QJ + 10Q^P^+3P'') I{3Qf+WQ^W+3P^) 

For ^ 1 in the expression for (S^) in Table I one obtains the same result as from Eq. (31) when substituting 
S" = 1/2 to 5/2 into it and replacing z by x. 

With = 1 + hy^, a general expression analogous to Eq. (31) is found 



^ [(25 + 1)P Qx] [Qx + Pf'+' + [(25 + l)P + Qx] (Q. - P?'+' 

^ ' 2F2[(Q^+p)2S+l_(Q^_p)2S+l] • ^ 
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Substituting S* = 1/2 to 5/2 into Eq.(34), wc obtain the expressions in the column 2 in Table I. Thus we suggest that 
Eq.(34) is the required correct expression for any S (we have no analytical proof for the general formula). 
We also suggest the following analytical expressions for the three correlations for any spin value 

US')') = (35) 
((S»)^) = ^S{S+l)hi, + (P^-3hl,)Q,S' ^ jj^j 



and 



It is seen that the exact relation 



{{Sn') = Iq^S-. (37) 



((5^)2) + {{syf) + {{S'f) = S{S + 1) (38) 

is satisfied. Also, when 5=1/2 to 5/2 is substituted into Eqs.(35, 36), we obtain the expressions in columns 3 and 4 
in Table I. 

Because in the present case the x and y-directions should be equivalent, exchanging x < — > y in Eqs. (34-37) and in 
Table I must also yield valid expressions. 

2.2. The general case 

We now take the 2;-direction as reference direction, supposing hz 0, and taking two equations from Eqs. (26) 

2(^2 — IhyzC-i = —iQzF-i^i — F-i,2 + , , 

Again, by specifying proper B operators, one obtains from these equations analytical expressions for the magnetization 
components and all correlations for each S. In Table II the explicit expressions for (5^), {{S^)^), {{S^)"^) and {{S^)^) 
for spin quantum number S = 1/2 to 5/2 are listed, where 

R^ = l + hlz + hi,. (40) 

Table II. The analytical expressions for (S"^), ((S'^)^), {{S'=f} and ((S'^)^) for S = 1/2 to 5/2. 

S (S^) ((5^)2) 

1/2 1/2Q^ 1/4 



1 

3/2 
2 

V2 2Q^ {3Qt+10Qi R^+3R*) 4(3Q| + 10Q^fl^+3fl«) 



4Q.(5Q^+3J?.^) WQt+eQl{R^+7)+6R' 
35QI+42QIr'^+3R'^ 35Q^ + 14Q^(3fl''+16)+3(fl''+32)fl^ 



S {{S^f) ((5^)2) 

1/2 1/4 1/4 

3QI + R:^ 3Ql+R^ 
bQl+R^ + 12hi^ 5Ql+R^ + 12hl^ 

I{Q^+W) 4(QfW) 

WQi+QQl{B^ + 7hl^)+liR^hl^ WQt+eQl{R' + 7hlJ+6hl^R' 



1 

3/2 
2 



r/r, 35Qt + 14Ql(3R^ + 16hlJ+3{R^+32hl^)R' 35Q^ + 14Qj (3fl" + 16/i^ J+3(fl"+32fe^ Jfl^ 
/ 4(3Q| + 10g2fl2_,_3^4) 4(3Q| + 10Q2_R2_|_3fl4) 
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and 



Similar to the case in subsection 2.1, we suggest that these quantities can be expressed by the following equations 
for any S. 

_ [(25 + 1)R - Q,]{Q, + + [(2^ + 1)R + Q,](Q, - R)^s+^ 

^(^.)2^ ^ 2S{S+1) + {R^~3)Q.S^ ^ ^^2) 
^^^2;-j2^ _ 25(5' + l)hl^ + (i?^ — 3hl.^)QzS^ ^^^^ 

((..)2. 2SiS + l)hl + {R'-3hl)Q.S^ ^ ^^^^ 

Here equation (38) is again satisfied. As long as (5^) is calculated, the other two components of the magnetization 
can be computed from the regularity conditions Eqs.(28). 

Because none of the three directions is special, we can also pick the x- or ^/-direction as reference direction. Therefore, 
in Eqs. (41-44) and in Table II, the formulas remain valid when performing the exchanges of coordinates, e.g., 
X ^ z ^ y ^ X. Note too, that the expressions of Eqs. (41-44) revert to those of Eqs. (34-37) of subsection 2.1 
after making the exchange x ^ y ^ z ^ x and setting hzx = 0. 

2.3. The effective (temperature-dependent) anisotropies 

The absolute value of the magnetization and the equilibrium angles and ipo of the magnetization direction are 
determined by the following equations: 

M^{T) = (5^)2 + {Syf + {S'f, (45) 



00 = arctan — , tpo = arctan — — . (46) 

Frobrich et alm^ have determined the effective (temperature-dependent) anisotropy coefficients K2ziT) non- 
perturbatively for the out-of-plane anisotropy case by minimizing the free energy with respect to the orientation 
angles of the magnetization. In this paper, we are also able to determine the coefficients K2x{T) and K2y{T) in the 
same manner. 

The free energy corresponding to the Heisenberg model with anisotropies in each direction reads 

F{T) = Fo{T) - K2x{T) sin^ e cos^ ip ~ K2y{T) sir? 9 s\T? ip - K2z{T) cos^ 9 

- B^M sm9 cos (fi - By Msm9smip - BzM cos9 (47) 

We have mentioned that, due to the relation (5^)^ -I- (S^)"^ + (5^)^ — 5(5 + 1), the three anisotropies are not 
independent of each other. Therefore we are allowed to transform for instance K2x K2x — ^2y, K2y 0, and 

K2z'^K2z~K2y. 

Choosing the xz-plane as reorientation plane by taking By = 0, which leads to (p — 0, the variation of dF/d9 = 
yields 

, BxM COS 9n — BzM sm9Q , , 

K2Z (T) - K2X (T) = °- . 48 

2 sm 00 cos 9o 

The result of the uniaxial out-of plane case— is obtained by taking K2x = 0. The angle 6*0 is the equilibrium angle. 

When considering for example the in-plane magnetization with K2x ^ and K2y ^ and K2z = (no out-of-plane 
anisotropy), one has two equations, 

dF{T)/d9^0, (49) 
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and 

dF{T)/dip:^0 (50) 

to determine K2x{T) and K2y{T) with the result: 

K2x[T) — — — , 

2 cos t^o 2 sm t^o cos 9?o 

K (T) BJ^__ByM__ 

' ~ 2 cos 00 2 sin 00 sin ^o' ' 



where the angles ^o, ^po are to be considered as equilibrium angles. 



III. NUMERICAL RESULTS 



In the numerical calculations, the strength of the exchange energy is set to J = 100. The parameters are scaled as 
in RefsilSiii, 

J ^ J/SiS + 1),B^ B/S, K2 K2/S{S - 1/2). (52) 

Scaling is introduced to obtain magnetization curves which arc approximately universal for different spin values S. 
The scaling of K2 guarantees the property \imT-ioK2{T) / K2 — 1. For the 3-dimensional (3D) case, we use a simple 
cubic lattice and for the 2-dimensional (2D) case, a square lattice. 



3.1. The 3D case 



For the 3D examples, we always put K2Z — and K2X — 1- For the uniaxial in-plane case {K2y — 0) the 
magnetization points in the x-direction. In the absence of magnetic fields, the corresponding scaled Curie point 
Tc/Tc{S = 1), normalized to that of 5 = 1, versus the spin quantum number S is plotted as squares in Fig.l. The 
scaled Tc slightly decreases with increasing S. 



CO 
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0.98 



0.96 



0.94 
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Wang etal. Fig.1 
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FIG. 1: The scaled Curie temperature Tc versus the spin quantum number S when there is only a single-ion anisotropy in one 
direction. Squares: 3D case. Upward triangles: 2D case without dipole-dipole interaction. Circles and downward triangles : 
2D case with dipole-dipole interaction (with strength g = 0.018), when the easy-axis is in plane and out of plane, respectively. 



We show some examples when including a magnetic field B = {Bx, By, B^) in addition to the anisotropics K2X = 
l,K2y ^ 0,K2z = 0. The direction of the magnetization will idtimately turn towards the field direction when 
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Wang et al. Fig. 2 
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FIG. 2: The angles of the direction of the magnetization (^o, V^o), defined by Eas. l)46^ . as functions of the temperature for 
K2z = 0, K2x = 1 and various K2y anisotropies and magnetic fields. The symbols are just to label curves. 



Wang et al. Fig .3 

1 90 




Temperature 

FIG. 3: The effective anisotropy K2xy{T) / K2xy{T = 0) and the equilibrium polar angle 6o for S — 1, 2 and 10 in the 3D case 
for in-plane anisotropy. The parameters are K2xy = K2x ~ K2y = 1 and = 0.1. 



the temperature goes above the Curie point (defined from {S^) = for the field-free case). For instance, when 
B = (0.1,0.1,0.1), (00, Vo) approach (57.74, 45) degrees, which is the [111] direction, and when B = (0, 0.1, 0), 
ifiQ approaches 90 degrees, as shown in Fig. 2. However, at low temperature, where the field is not strong enough 
to overcome the effect of the anisotropy, the magnetization will be close to the direction along which the anisotropy 
is strongest. With increasing temperature, the magnetization gradually turns towards the field direction. At zero 
temperature, the initial angles of ipo are 3.60 and 2.74 degrees for K2y ~ 0.5 and 0, respectively, for both apphed 
fields B = (0.1, 0.1, 0.1) and (0, 0.1, 0). When K2y = 0.5, the magnetization turns more easily than with K2y ~ 0. 
Note that the dashed and dot-dashed lines are very similar. This reflects the fact that the magnetization component 
in the plane is almost the same for the two cases of K2y = 0.5 and 0. If B = (0.1, 0.1, 0), the magnetization is in the 
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FIG. 4: The components of the magnetization (a) {S'')/S and (b) {S^) / S as functions of the external field in z-direction Bz 
for S = 1, 2 and 10 and at temperatures T = O.STc and O.STc. Tc is the scaled Curie temperature of the 3D case shown by 
squares in Fig.l. 



xy-plane (^o = 90"). The curve of the angle ipo versus temperature is identical to the solid line in Fig. 2 for K2y = 0.5 
and to the dotted line for K2y — 0. 

As a further example we discuss the in-plane anisotropy case, with if 2x3/ = -?^2x = = 1. The effective anisotropy 
coefficient as a function of the temperature can be calculated by Eq. (51). Figure 3 shows examples for the spin 
values S = 1, 2 and 10 calculated for = 0.1 and = By = 0. At zero temperature , K2xy{T = 0) is the same as 
the anisotropy parameter of the Hamiltonian K2xy With increasing temperature, K2xy{T) decreases. Close to the 
Curie point, the coefficients approach zero. With the present choice of parameters the reorientation temperature is 
close to the Curie point; in general they can be quite different. In the intermediate temperature range, one observes a 
sensitivity to the spin value: the larger the spin S, the smaller the effective coefficient. The larger the spin the more 
increases the range of a linear behaviour of the effective anisotropy, which is in accordance with mean field results 
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where for 5 — > oo or for classical spins one obtains a linear behaviour K2{T) = K2{0){l — T /Tc)^. Figure 3 also shows 
the corresponding equilibrium angle defined by Eq.(46). The in-plane anisotropy should force the magnetization 
into the xy- plane. However, because of a small field in z-direction Bz = 0.1, the magnetization deviates slightly 
from the plane at zero temperature, and there is a small initial angle. It is the same for each S. With increasing 
temperature, the magnetization reorients slowly, owing to the action of the field. Close to the Curie point, the effective 
anisotropy approaches zero, so that the magnetization reorients rapidly in the direction of the z-axis, the direction of 
the external field. 

The orientation angle Oq of the magnetization depends on the temperature as well as on the external field. In 
Figure 4, we plot the magnetization components {S'^)/S and (5^}/S' versus the unsealed external field at two 
temperatures for S = 1, 2 and 10. Because of scaling the temperature dependence for different spin values looks 
similar. The reorientation fields at which the reorientation to the field direction has occurred however are different 
and are given by Bz = 0.76, 0.38 and 0.08, respectively at T = 0.8Tc. The larger the spin value 5', the easier is 
the reorientation of the magnetization. To compare the reorientation behaviour at different temperatures, we also 
show the results at T = 0.5Tc. In this case, the magnitudes of the field at which the magnetization reorientation 
occurs are Bz = 0.98, 0.49 and 0.10, respectively, which are higher than those at temperature T = O.STc. At a lower 
temperature, it is more difficult for the magnetization to turn under the action of the field, as one expects. 



3.2. The 2D-case (monolayer) 

For the discussion of a monolayer (the 2D-case), we add the dipolc-dipole interaction (DI) to the Hamiltonian of 
Eq.(l) 




where is the distance between sites i and j, and its components along three axes are Tijufj, '>~ij''J'ij and rijuf^, 
respectively. The higher-order Green's functions containing DI terms are also decoupled by RPA, i.e. by Eq.(lO). 
As a result, there appear additive terms due to DI in of Eq.(16), which include non-dispersive parts as well as 
dispersive parts. As an approximation, we neglect the dispersive parts and retain the non-dispersive parts only. Such 
an approximation was used before in Ref.^"^. The full RPA treatment of the dipole coupling would require complex 
arithmetic numerically (see Appendix A of Refii^). An exception is the uniaxial out-of-plane case, for which it was 
shown in this Appendix that the mean-field results are close to RPA. Using complex arithmetic numerically is not in 
principle a hindrance, but is quite tedious to perform. One obtains 

He. = (5") (Jo - Jk) +B^+ i^2a(^")*a - .9(^") J] ^[1 ^ 3^^], « = I, J/, z. (54) 

In the last term of Eq.(54), the subscript means the origin, and the summation covers all lattice sites except for the 
origin. The lattice sum can be done and Eq.(54) becomes explicitly 

H^,y = {S'''y){Jo-Ju) + B,.,y + K2^,y{S^^y)<f^,y + gTo{S^'y), 
Hz = {S'){Jo-A)+Bz + K2z{S')<fz-2gTo{S'), (55) 

where Tq — J2i m (;2+to2-)5/2 = 4.5165. To compare with the results in Ref;— , we choose g — 0.018 (a value which 
corresponds to Ni with the present set of parameters) in most cases. Note that the DI strength g is also scaled as 



g^g/S{S+l). (56) 

Correspondingly, the expressions for the effective anisotropics Eqs.(48, 51) will have additional terms coming from 
the DI, as in Refpi^ for the out-of-plane case. 

Now we proceed to some numerical examples. Owing to the Mermin- Wagner theorem.^'*, a Heisenberg monolayer 
without some kind of anisotropy has no spontaneous magnetization. A spontaneous magnetization can be achieved 
by DI alone. Figure 5 shows the absolute value of the in-plane magnetization due to the action of DI as function of 
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the temperature. One observes that for different spin quantum numbers S, the scaled Curie temperature Tc is the 
same for a fixed DI strength. The inner panel shows the T^ — g curve in the range oi g — 0.01 — 1. When g approaches 
zero, Tc decreases drastically, and approaches zero for g ^ in accordance with the Mermin- Wagner theorem. 



Wang et aL Fig.5 




Temperature 



FIG. 5: Magnetization-temperature curve for 5 = 1,2 and 5 for the monolayer with DI strength g — 0.018 and 1, respectively. 
M is the magnetization defined in Ea. H45^ for the in-plane situation. When g is fixed, Tc is same for any S. The inner panel: 
Tc as a function of the DI strength g. 

If there is only a single-ion anisotropy in addition to the Heisenberg exchange energy term, a spontaneous magne- 
tization of the monolayer is also induced. The upward triangles in Fig.l show the scaled Curie temperature Tc as a 
function of S when there is an in-plane anisotropy in one direction. The Curie temperatures are lower and decrease 
more rapidly than those for the 3D case. The Curie points for the 2D case are lower than those for the 3D case for 
the same S, since the number of nearest neighbours of each site in a 2D square lattice is four, while that in a cubic 
3D lattice is six. 

With an external field in addition to the single-ion anisotropies, the magnetization will rotate. The rotation in the 
monolayer case exhibits (except for the temperature scale) the same features as those of the 3D system. If we plot 
the 6q, ifQ — T curves, the qualitative behaviour will be the same as in Fig. 2. With the same parameters, the initial 
angles are the same as in Fig. 2. If, for example, B = (0.1, 0.1, 0.1), K2z = and = 1, the angles {Oq, (po) at zero 
temperature are (87.26, 3.61) degrees for K2y = 0.5 and (87.26, 2.64) degrees for K2y = 0. Above the Curie point, 
{9o, (fo) again approach the field direction (57.74, 45) degrees, i.e. the [111] direction. 

Including both the single-ion anisotropy and DI, we plot in Fig. 1 the scaled Curie temperature Tc as a function of 
S for the uniaxial cases K2x = 1 (in-plane) and K2z = 1 (out-of-planc) using g — 0.018. It is seen that, if the easy-axis 
is along the out-of-plane direction, the effect of DI is to lower Tc- If the easy-axis is in the in-plane direction, the 
effect of DI is to raise Tc. This is clear from Eqs.(55) where the effect of DI can be understood as a renormalization 
of the non-dispersive exchange interaction terms: 

Ms--'>^) ^ iJo + gTo){S-'-y), 

MS') ^ iJo-2gTo){S'). (57) 

Jk in Eq.(55) is kept unchanged. As a result, the exchange energy strength along the z-direction is suppressed and 
that in the xy-plane is enhanced. The figure shows that in the out-of-plane case, the Curie temperature Tc decreases 
(compared to g = 0) more strongly than it increases for the in-plane case. 

Now we investigate the in-plane (or easy-plane) anisotropy case, where K2xy = K2x = K2y — 1. Figure 6 shows 
the effective anisotropy coefficient, defined by Eq.(51), as a function of the temperature for = 1, 2 and 10, when 
Bz = 0.1 is taken. At zero temperature, K2xy{Q) is the same as the parameter K2xy The effective anisotropy decreases 
with increasing temperature, and at the Curie point it approaches zero. In the intermediate temperature range it is 
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Wang et al. Fig.6 
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FIG. 6: Effective anisotropy K2xy{T)/ K2xy and equilibrium polar angle 60 for S = 1, 2 and 10 in the 2D case. The parameters 
are same as in Fig. 3. 



the smaller the larger the spin quantum number S is. Figure 6 also shows the equilibrium polar angle defined by 
Eq.(46) when — 0.1. All the properties in Fig.6 are similar to those in Fig. 3. The only difference is that the Curie 
temperatures in the 2D case are by more than a factor of two smaller (Fig.l) and more different for the various spin 
values (Fig.6) than in the 3D case (Fig. 3). 

The polar angle of the magnetization depends on the temperature as well as on the external field. It is unnecessary 
to plot the components of the magnetization as functions of the external field at various temperatures as in Fig. 4 
(the 3D case) because they are qualitatively so similar. 

In RefsfiSiii, the case of the out-of-plane anisotropy was studied in detail. Since the method developed in the 
present paper also enables the study of the in-plane anisotropy, it is natural to compare both cases. We mention 
the following differences, which can all be explained on the basis of Eqs. (55). The first concerns the case where 
there is no external field. For an out-of-plane anisotropy, there is a magnetization reorientation from the z-direction 
to the xy-plane due to DI even without an external fielcU^. At zero temperature, the magnetization points to the 
^-direction because of the out-of-plane anisotropy term K2z- With increasing temperature, the magnetization turns 
to the xy-plane owing to the DI. In the case of in-plane anisotropy, there is no spin reorientation in the absence of 
the field, because the magnetization is already in the xy-plane at zero temperature. 

Secondly, with the same DI strength and the same field (perpendicular to the anisotropy direction), the magneti- 
zation reorientation temperature is different for the two anisotropy cases as shown in Figure 7. In Fig. 7(a) and (b), 
the solid lines are for 5=1 and the dotted lines are for S — 10, both with g — 0.018. Figure 7(a) depicts the in-plane 
anisotropy case with an external field = 0.1 along the z-direction, and figure 7(b) the out-of-plane anisotropy case 
with an external field B^ =0.1 along the cc-direction. The magnetization reorientation temperature of the latter case 
is lower than that of the former case. This can be understood from Eqs. (55), from which one sees that the dipole 
interaction enhances the effect of the in-plane anisotropy term, whereas the opposite is true for the out-of-plane case; 
i.e. the role of DI can be considered as a renormalization of the anisotropy terms in the following sense 

if 2.$. ^ K2,^,-2gT^. (58) 

This also explains, thirdly, the effect of the variation in the strength of the DI, as shown in Figure 7(a) and (b). In 
the case of in-plane anisotropy, increasing the DI strength hinders the reorientation from in-plane to out-of-plane, see 
Fig. 7(a), whereas the reorientation from out-of-plane to in-plane is enforced in the case of the out-of-plane anisotropy. 
It is seen that the DI suppresses the out-of-plane anisotropy and enhances the in-plane anisotropy. 
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Wang et al. Fig. 7 




Temperature 

FIG. 7: (a) Effective anisotropy K2xy{T)/K2xy and equilibrium polar angle 60 for 5* = 1 and 10 for the monolayer with DI 
strengths g — 0.018 and 0.15 for the in-plane case. The parameters are K2xy = 1 and — 0.1. (b) Effective anisotropy 
K2z{T)/K2z and equilibrium polar angle 9o for 5* = 1 and 10 for the monolayer with DI strengths g = 0.018 and 0.08 for the 
out-of-plane case. The parameters are K2z ~ 1 and Bx ~ 0.1. 



The fourth difference involves the variation of the external field. If an external field is applied along the a;- 
direction in the case of an out-of-plane anisotropy, the magnetization reorients easier with DI than without DI. This 
is shown by comparison of the dashed and solid lines in Fig. 8. The reason is that we can also interpret the effect of 
the DI as a renormalization of the external field, as has been done in Ref>i^, see Eq. H55|) 

^ -2.gTo(5"). (59) 

This means that the effect of the DI is to enhance a magnetic field in the xy-plane, and to weaken an external field 
along the z-direction. Therefore, if an external field Bz is applied along the z-direction in the case of an in-plane 
anisotropy, it is more difficult for the magnetization to reorient into the out-of-plane direction with DI than without 
DI. Comparison of the dotted and solid lines in Fig. 8 provides an example. 

IV. CONCLUSIONS AND DISCUSSION 

In this paper, we have treated a Heisenberg Hamiltonian with single-ion anisotropics in all three directions of space 
by applying many-body Green's function theory, thus generalizing previous work dealing with uniaxial out-of-plane 
anisotropy. To do so, we selected the spin operators 5*^, 5*^ and as basis operators to construct the Green's functions 
instead of the commonly used operators 5*+, S~ and . We stress that in the derivation of the formalism none 
of the three axes is special, so that we are always able to select a reference direction along which the magnetization 
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FIG. 8: The components of the magnetization (S^) and {S^) as functions of the external field for 5 = 1 at the temperature 
T/Tc = 0.8. The dashed lines correspond to the out-of-plane anisotropy with a field in a;-direction. The solid and dotted 
lines correspond to the in-plane anisotropy with a field in z-direction. If K2xy = , K2z = 1 and g = and a field B^, the 
corresponding magnetization curves are identical to the solid lines with (S^) and (S^) exchanged. 



component is not zero. The key is the generalization of the Anderson-Callen decoupling to the terms from the single- 
ion anisotropy in all directions of space, Eqs.(ll, 12). A formal advantage of the use of the present operators is that 
the matrix governing the equations of motion is Hermitian, whereas the ususal formalism leads to real non-Hermitian 
matrices, which are more difficult to treat. A main result of the paper is that analytical expressions can be derived for 
the three components of the magnetization and the correlations ((S*^)^), {{S'^)^) and {{S^)^) for any spin quantum 
number. The formalism obtained in this paper enables a treatment of both the in-plane and out-of-plane anisotropics. 
Numerical calculations were carried out for various parameters for both 3D and 2D systems. 

For the 2D system, we included the dipole-dipole interaction. Its effect on the magnetization can be explained on 
the basis of Ea. (|55|l . and can be interpreted as a renormalization of the nondispersive part of the exchange interaction 
terms, or of the single- ion anisotropy terms or of the external field, see Eas. (|57l 1581 159() . The reorientation of the 
magnetization from in-plane to out-of-plane under an external field occurs more easily without the DI than with it, 
or in other words the DI supports the reorientation from the out-of plane to the in-plane direction. 

Because of Eq.(55), may appear negative. In this case the regularity condition Eq.(28) would be meaningless 
because the z-component magnetization should be within < (S^) < S. Only when the value of g is under a 
certain value can there is a non-zero (S^) without a field in z-direction. From Eq.(58), the condition should be 
gTo < S{K2z^z - 2gTo). Then one should have 2SK2z (l - i^[S[S + 1) - (S^S*^)]) > gTa{l + 25). Since (5^5^) > 
S{S + l)/3, one has g < 2K2z{2S - l)/3To(l + 2S) = 0.1476(25 - 1)/(1 + 25), where = 1 is taken. The upper 
limit increases with spin quantum number. However, because we use scaled parameters in this paper, Eq.(52), we 
have g < 5(5 1)2^2^(25 - l)/3To(l-f 25)5(5 - 1/2) = (5 + l)4/3ro(l + 25) = 0.2952(5 -f l)/(25-H 1). The upper 
limit decreases with spin quantum number. When the value of g is over this limit, the z-component of magnetization 
has to be zero. Figure 5 provides an example where the magnetization lies in the plane. 

Finally, we point out that the numerical results are all calculated with scaled parameters in order to obtain universal 
results, in accordance with Refs>i2ii^. With increasing 5, the scaled exchange J and anisotropy K2 become smaller 
as can be seen in Ea. (|52|l . This results in the conclusions in Section 3 that, with increasing 5, the following quantities 
decrease: the scaled Tc in Fig.l, the effective anisotropy coefficients in Figs. 3, 6 and 7, and the field under which the 
reorientation occurs in Fig. 4. If the parameters are not scaled but kept unchanged, the conclusions are the other way 
round. A larger 5 means a larger magnetization. It will, therefore, result in a higher Curie temperature. Consequently, 
the effective anisotropy coefficients are larger for larger 5 at a given temperature. For instance, at zero temperature, 
the effective anisotropy coefficient should be K2S{S — 1/2). A larger magnetization is more difficult to turn by an 
external field. As a result, the magnitude of the field at which the magnetization reorientation occurs should rise with 
increasing 5 for a given temperature. 
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